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Universal collision rate constants are calculated for ultracold collisions of two like bosonic or fermionic heteronuclear alkali- 
metal dimers involving the species Li, Na, K, Rb, or Cs. Universal collisions are those for which the short range probability of 
a reactive or quenching collision is unity such that a collision removes a pair of molecules from the sample. In this case, the 
collision rates are determined by universal quantum dynamics at very long range compared to the chemical bond length. We 
calculate the universal rate constants for reaction of the reactive dimers in their ground vibrational state v = and for vibrational 
quenching of non-reactive dimers with v > 1 . Using the known dipole moments and estimated van der Waals coefficients of 
each species, we calculate electric field dependent loss rate constants for collisions of molecules tightly confined to quasi-two- 
dimensional geometry by a one-dimensional optical lattice. A simple scaling relation of the quasi-two-dimensional loss rate 
constants with dipole strength, trap frequency and collision energy is given for like bosons or like fermions. It should be possible 
to stabilize ultracold dimers of any of these species against destructive collisions by confining them in a lattice and orienting 
them with electric field of less than 20 kV/cm. 
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1 Introduction 

Quite spectacular success has been achieved in recent years 
in working with gases or lattices of ultracold atoms cooled to 
temperatures on the order of 1 fiK or less. This has permitted 
the achievement of quantum degeneracy with either bosonic or 
fermionic isotopes of various atomic species where the ther- 
mal de Broglie wavelength becomes of the same order as or 
larger than the mean distance between atoms. Quite precise 
control of the various properties of such systems is possible 
through state selection, trap design, and magnetic tuning of 
scattering resonances. A number of reviews or books cover- 
ing this work have appeared, including topics such as Bose- 
Einstein condensation, 1 2 the quantum properties of fermionic 
gasesj2El and magnetically tunable Feshbach resonance con- 
trol of collision properties.SEl j n addition, much work has 
been directed towards creating lattice structures of such cold 
atoms using standing wave light patterns to make optical lat- 
tices of various geometric configurations 

Recent work has now succeeded in making ultracold 
molecules in their stable vibrational and rotational ground 
electronic state with temperature T < 1 /J.K by using magnetic 
and electromagnetic field control of ensembles of ultracold 
atoms. This has been achieved for 40 K 87 Rb fermions^ and 
133 Cs2 bosons 15 , and for 87 Rb2 bosons in the collisionally un- 
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stable v = level of the lowest 3 E+ stated These successful 
experiments used a two-step process to make the molecules. 
First, magneto-association of two very cold atoms makes a 
very weakly bound "Feshbach molecule", 8 which is then con- 
verted by a coherent STIRAP process to make a ground state 
molecule in its v = 0, / = ground state in a single state of nu- 
clear spin, where v and J are the respective quantum numbers 
for vibration and rotation. This builds on pioneering earlier 
work in much more dilute cold atomic gases around 100 /J.K to 
make ground state RbCsES Cs2,^or LiCs^ molecules. Cold 
molecules open up many new opportunities for study,^2H2] 
since they have more complex internal structure and different 
long range potentials. If dipolar, their properties and collisions 
can be controlled by an electric field in addition to magnetic 
or electro-magnetic fields J 23 l 24 l 

Having ultracold molecules also introduces the possibility 
of chemistry and reactions at ultralow energies, with precise 
control of the initial internal states and translational energy 
of the reactantsJ 22 ! 25 ! Having such collisions is good if one 
wishes to study such chemistry following the precise prepa- 
ration of the initial states. 26 Collisions are bad, however, if 
one wishes to keep the molecules for simulating condensed 
matter systems or doing complex control like quantum com- 
puting, since reactive collisions can rapidly remove molecules 
from the gas in question. Dipolar molecules offer some spe- 
cial features such as the possibility of orienting them by an 
electric field. Reaction rates can be strongly suppressed if the 
molecules are oriented to have repulsive dipolar interactions 
while confined to move in a two-dimensional (2D) plane by 
an optical lattice wave guide. Such suppression of reaction 
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rates in this quasi-2D geometry has been both predictecPZH^H 
and demonstrated for 40 K 87 Rb. 30 Controlling internal spin can 
also be used to decrease reaction rates in the case of fermionic 
molecules like 40 K 87 Rb, since identical fermions can only col- 
lide via odd partial waves that have centrifugal barriers to re- 
action.EMl 

Here we examine some basic aspects of ultracold chem- 
istry of highly reactive molecules. We restrict our consider- 
ations to the special case that the two interacting molecules 
have a unit probability of a chemical reaction or an inelastic 
quenching collision if they approach one another within typi- 
cal chemical interaction distances, on the order of 1 nm or less. 
Thus, we are considering the quantum threshold limit to the 
standard Langevin model. 32 34 In this highly reactive limit the 
long range potential between the molecules controls how they 
get together subject to experimental control on an ultralow en- 
ergy scale. This paper will examine various universal aspects 
of reaction rates that are solely gover ned by such long range 
interactions of two reactive moleculesJ 28 ! 34 ! Such interactions 
are sensitive to the Bose/Fermi character of the molecules, to 
the traps used to confine them, and to electric, magnetic, and 
electromagnetic fields in the case of polar molecules. We will 
consider, in particular, the ten different molecular dimers com- 
prised of two alkali metal atoms from the group Li, Na, K, Rb, 
and Cs. 

The theory will be confined to the lowest temperatures 
where only the lowest partial waves allowed by symmetry con- 
tribute. We will apply both analytic and numerical approaches 
to molecules colliding in the 3D geometry of free space or in 
the quasi-2D geometry of a confining optical lattice. The basic 
theory related to ultracold molecules, their formation, states, 
collisions, dipolar properties, and response to external fields is 
discussed in detail in the introductory book by Krems et al. 35 
and the Faraday Discussions 142.^ We will review this the- 
ory as needed here and apply it to dipolar mixed alkali dimer 
molecules. 



vides an efficient and effective tool to make near threshold 
bound states known as Feshbach molecules, 8 9 with binding 
energies E/h less than 1 MHz. These can then be optically 
converted via a coherent Raman process to much more d eeply 
bound vibrational levels,E2El including the ground stateM 4 -^ 
Typical temperatures are well below 1 jiK and densities can 
be on the order of 10 12 molecules/cm 3 or even larger. Since 
this can be done by coherent quantum dynamics that causes no 
heating, the molecules have the same temperature as the ini- 
tial atoms. It is even possible to make an optical lattice array 
of many single trapping cellsj 15 l 41 lin which exactly two atoms 
are present. Upon associating the atoms, one then has an array 
of molecules, each of which is trapped in its own lattice cell, 
with confinement possible to tens of nm and intercell spacing 
of hundreds of nm. 

In principle any diatomic molecule could be made that is 
comprised of atoms that can be cooled and trapped, including 
species such as Li, Na, K, Rb, Cs, Ca, Sr, Yb, or Cr. In practice 
so far, ultracold molecule formation has been concentrated on 
alkali-metal atom dimers. Among the alkali-metal-species Li, 
Na, K, Rb, and Cs, it is known 42 that five mixed dimers have 
exoergic reactive collision channels when in their vibrational 
and rotational ground state, v = 0, J = 0, namely, LiNa, LiK, 
LiRb, LiCs, and KRb. The other five, NaK, NaRb, NaCs, KCs, 
and RbCs, have no reactive channels for v = 0, J = 0. 

A molecule in a pure rotational eigenstate has no dipole mo- 
ment. Calculating the energy of a molecular dipole in an elec- 
tric field F is explained in Chapter 2 of the book by Krems 
et al. We apply this method to the case of a 'E + sym- 
metric top rotor in vibrational level v in its ground rotational 
state J = M = 0, where M is the projection of total angular 
momentum J along the direction of F. Upon expanding the 
wave function in symmetric top basis states \JMAv) = |/00v), 
where the body-frame projection A = 0, the energy E g {F) of 
the lowest energy ground state \g(F )) of the manifold is found 
by diagonalizing the Hamiltonian matrix H moi with matrix el- 
ements: 



2 Collisions in free space 
2.1 Ultracold polar molecules 

While general cooling methods such as buffer gas cooling and 
Stark deceleration to load a molecu le trap are being developed 
for a variety of moleculesj 21 ! 22 ! 35 ! such schemes so far have 
been restricted to temperatures significantly above 1 mK and 
very low phase space density, many orders of magnitude re- 
moved from quantum degeneracy. To date, high phase space 
density has been restricted to molecules that can be made di- 
rectly by associating two atoms that are at or near the quan- 
tum degenerate regime of temperature and density. While 
initial proposals to do this involved photoassociation of the 
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where B v = /z 2 (000v| 1 /r 2 |000v) /M is the rotational constant 
for the molecule with mass M, r is the interatomic separation, 
and d m is the body-frame permanent molecular dipole mo- 
ment. The energy E g (F) approaches the energy of the v,J = 
molecular level as F —> 0. The field-dependent dipole moment 
for the field-dressed ground state is d(F) = —dE g (F)/dF, 
which approaches d(F) — > d m (Fd m )/(3B V ) as F — » 0. 

In order to make estimates for the dipolar collision prop- 
erties of the ground and lowest vibrational levels of the ten 
mixed-alkali-dimer species, we will take the dipole moments 
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d m = d e and rotational constants B e for the ground E + state 
evaluated at the equilibrium internuclear distance R e from the 
calculations of Aymar and Dulieu.ElFigure[I]shows the field- 
dependent dipole moments d(F) as a function of F for the 
respective reactive and non-reactive sets of species. 

A key aspect to note from Fig.[TJis that for all species except 
LiNa and KRb, dipole moments on the order of 0.4 au (1 D = 
3.336 x 10~ 30 Coulomb meter) can be achieved at relatively 
modest electric fields of 10 kV/cm or less. As F increases, 
d(F) approaches and eventually reaches the magnitude of the 
molecular dipole d m . This is not yet achieved at the 20 kV/cm 
maximum in Fig. [T] The KRb experiments on collisions in 
an electric field were carried out with F up to 5 kV/cm and 
were able to reach dipole moments only on the order of 0.08 
au (0.2 D)J2ffil Except for LiNa, which has an even smaller 
dipole moment than KRb, all of the other species should be 
capable of being used with dipole moments on the order of 
0.4 au (ID) or more to enable quite interesting experiments to 
control threshold collision dynamics. 

2.2 Characteristic energy and length scales 

In order to get a good understanding of the properties of ul- 
tracold molecule interactions, it is important to understand 
the various length scales L and corresponding energy scales 
h 2 j(2\ih 2 ) associated with ultracold phenomena. These are 
the de Broglie wavelength 2n/k and the lengths R e , a, a/,, 
and cid associated with the respective chemical, van der Waals, 
harmonic, and dipolar interactions. For a relative collision ki- 
netic energy of E — h 2 k 2 / (2ju), where jx = M/2 is the reduced 
mass of the molecule pair of mass M, we define a characteris- 
tic energy-dependent length a^= 1/k. The characteristic van 
der Waals length is a = [27z;/r(l/4) 2 ] (Z/iQ/ft 2 ) 1 ' 4 , where 
— C(, /R 6 is the van der Waals dispersion potential, R is the in- 
termolecular separation, and 2n/T(l /4) 2 « 0.47799.13221 The 
harmonic length for trapping frequency £2 is the characteristic 
length of ground state motion a/, = Wh/ (juQ). The potential 
energy of interaction between two molecular dipoles separated 
by distance R is d(F) 2 (l — 3 cos 2 9)/R i , where 9 is the angle 
between dipole orientation and the intermolecular axis. The 
dipole length a c i(F) = jid{F) 2 /h 2 is defined to be the length 
where h 2 /(iia 2 ) = d(F) 2 ja\. 

In order to estimate a, we need the Cg value for the various 
species. For a polar molecule there are two contributions to the 
effective C(, interaction. One is the electronic contribution Cf, 
due to the second-order response through excited electronic 
states. A simple Unsold approximation gives Cf = (3/4)C/a 2 , 
where U ps 0.055(7) atomic units is a mean excitation energy 
and a is the dipole polarizability.^MaThis gives a magnitude 
on the order of Cf 10 4 au (1 au = E h a% = 9.573 x 10~ 26 
J nm 6 ) for the mixed alkali dimers. For most species, the 
much larger and dominant contribution to Cs is the rotational 
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Fig. 1 Dipole moment d(F)/D versus F for the five reactive mixed 
alkali-metal-species (upper panel) and for the five non-reactive 
mixed alkali-metal-species (lower panel), where D = 0.3934 au = 
3.336 x 1(T 30 Cm. 

dipole part with a magnitude Cg ot = d^/(6B v )&l Thus, for 
estimation purposes for the vibrational ground state, we use 
C(, ~ (3/4)Ua 2 + d^/ (6B e ). This approximation can be com- 
pared to the calculations of Kotochigova 47 for KRb and RbCs, 
where our estimates give respectively 10500 + 2800 au and 
15000+ 113000 au for the two contributing terms, giving C(, 
values 10 and 20 per cent less than the corresponding ab ini- 
tio calculations. Since the length a is not very sensitive to C^, 

1/4 

varying only as C 6 , this approximation makes small errors 
on the order of 5 per cent or less in our qualitative estimates 
of a. Consequently our estimates for universal van der Waals 
rate constants are likely to be in error by 15 per cent or less for 
all species except perhaps LiNa, which has the smallest dipole 
moment. 

Figure|2]illustrates the various characteristic lengths for the 
ten mixed alkali dimers. The shortest distance is the chemical 
bond length R e , which is much less than 1 nmj^with a corre- 
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Fig. 2 Characteristic length scales R e (black dots), a (squares), a/, at 
Q. = 27r(30kHz) (blue dots), thermal expectation value of l/k at 200 
nK (crosses), and a c [ (stars for d = 0.4 au (1 D) and diamonds for 
d = d m ) for the ten mixed alkali-metal species. 



sponding energy scale E/h on the order of the chemical bond 
energy of 100 THz. The next largest energy scale, on the order 
of MHz, is that of the van der Waals length a, which ranges 
between 6 nm for KRb and 30 nm for LiCs. The harmonic 
confinement by an optical lattice typically has H/(2n) on the 
order of tens of kHz. The Figure shows that the confinement 
ah is on the order of 100 nm for Q. = 27t(30kHz). Figure [2] 
also shows the thermal average of 1 jk is on the order of a few 
hundred nm and is larger than either a or a/, when collision 
energy T is 200 nK, or k B T /h = 4 kHz, a typical value for 
ultracold systems. Note that the product ka^ = (2E / (hGl)) 1 ' 2 
is independent of the mass of the species for a given trap fre- 
quency Q.. The trap-induced confinement is less than 1 jk 
when E < \t£l, or E/k B < 700 nK, for an H = 27Z;(30kHz) 
trap. 

The dipole length depends on the electric field F. Fig- 
ure|5]shows a^F) for two values of the dipole strength d(F), 
namely, for 0.4 au (1 D) if this magnitude is possible, and for 
its maximum allowed value d m . These dipole length scales, 
respectively on the order of 1000 nm and 10000 nm, are the 
largest length scales in the problem, with the exception of 
LiNa, which has the smallest dipole moment. At d(F) = 0.4 
au, the corresponding energy scales are on the order of 100 
Hz, but over 1 kHz for LiNa, LiK and NaK. 

2.3 Universal collision rate constants 

The basic theory of ultracold collisions is covered in Chap- 
ters 1 and 3 of Rrems et al. 35 and by HutsonEMl and Chin et 
alM. We will summarize here the essential points using the 
formulation with quantum defect theory (QDT) by Idziaszek 
and Julienne^ElFor the moment we will ignore the role of in- 



ternal spin structure in the ground molecular stateJSEI] since 
we are considering highly reactive collisions where this seems 
not to be relevant. 31 The essence of the QDT approach is the 
separation of energy and length scales, as seen in Fig. [2] The 
overall effect of the short range "chemistry" zone with R <C a 
is summarized in the outer van der Waals zone by two dimen- 
sionless QDT parameters s and < y < 1, which respectively 
represent a short range phase and the chemical reactivity in 
the entrance channel j; the scattering length of the entrance 
channel is parameterized by s, which can take on any value, 
and y = 1 represents unit probability of short range loss from 
the entrance channel, whereas y = means no loss. In the 
ultracold domain, we only need to consider the lowest partial 
wave j to index the channel, where j = for s-waves for like 
bosons or for unlike bosons or fermions, and j = 1 for p-waves 
for like fermions or for unlike bosons or fermions. 

The case of y = 1 corresponds to a special "universal" case, 
where there is unit probability of short range loss. In this case, 
the elastic and inelastic or reactive collision rates depend only 
on the quantum scattering by the long range potential at dis- 
tances R > a. Any incoming scattering flux that penetrates 
inside a experiences no reflection back into the entrance chan- 
nel since it is lost to reactive or inelastic channels. No scat- 
tering resonances can exist in this case. The rate constant for 
loss from the entrance channel Kj{E) for a relative collision 
kinetic energy E for a van der Waals potential takes on the 
following very simple universal form, independent of the de- 
tails of the short range potential or dynamics, at low collision 
energy where kd <C 1 , 
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where the identical particle factor g = 2 if the collision part- 
ners are identical bosons or identical fermions, and d\ = 
ar(i) 6 /(1447T 2 r(^) 2 ) 1.064aJSIlf the colliding particles 
are not identical, then the rate constant is K ls = Kq +K\ s (E) 
with g = 1, and the second term becomes much smaller than 
the first as k — > 0. According to the thresold law for a van der 
Waals potential, the thermal average is independent of tem- 
perature T for the s-wave but varies linearly with T for the 
/?-wave, 



Kf(T) 



where the factor 

• i ^^6 /r/3\2 



\5l3d s k B T/h, (4) 

for identical fermions comes from 
r(5) 6 /r(|) 2 ~ 1513. The density n of a uniform gas of iden- 
tical bosons or fermions varies as h = -Kf(T)n 2 . In the 
case of a two-species gas, the two densities ti\ and «2 vary 
as h\ = ii2 = —K^n\ni- 

The universal rate will only apply to species with y = 1 and 
thus to all five of the highly reactive species in the upper panel 
of Fig. [T] This has been demonstrated for KRb for T < 1 
jUK^ However, there is good evidence that excited vibra- 
tional levels of alkali dimer molecules collisionally quench to 
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Fig. 3 Universal rate constants for like bosons, A^ s ( red diamonds), 
like fermions K\ s (blue dots), and the unitarity limit, Kj u (green star 
for j = 0, green plus for j = 1), for van der Waals collisions at 
T = 200 nK for the ten mixed alkali-metal species. These are 
expected to give the correct order of magnitude of for reaction of 
v = 0,7 = for the five reactive species LiNa, LiK, LiRb, LICs, 
KRb and for vibrational quenching of low v > 1 levels for the five 
species that are not reactive for v = 0, NaK, NaRb, NaCs, KCs, and 
RbCs. As r^o, A-[ s (r) varies linearly with T, Kf u (T) varies as 

1 /y/T, and Kq is independent of T. As T increases the universal 
rate limits only apply as long as K l ?(T) < K^ U (T). 



Consequently the universal rate constants only apply if they 
are smaller than this limit, requiring ka < 1/4 for j = and 
fc 3 a 3 < 1/4 for j = I. The corresponding p-wave unitarity 
limit K\ U (T) for like fermions is 3 times larger than Kq u (T) 
for like bosons. 

All five alkali-metal species have bosonic isotopes, but only 
Li and K have stable fermionic isotopes. Consequently sta- 
ble NaRb, NaCs, and RbCs fermions do not exist, but 22 Na 
and 134 Cs have half lives longer than 2 years. The fermionic 
molecules tend to have rate constants much less than the uni- 
tarity limit at this T. The bosonic molecules are closer to their 
unitarity limit for 200 nK, especially those with the highest 
dipole moments. 

Since the decay rate per molecule for a gas of density n 
is K^n, the lifetime of the molecule is x — (Kpi) -1 . Thus, 
taking n = 10 12 cm~ 3 as a typical ultracold gas density, one 
gets lifetimes of 1 s and 1 ms for respective rate constants 
of 10~ 12 cm 3 /s and 10 -9 cm 3 /s. Thus reactive bosons and 
unlike fermions will have typical lifetimes of a few ms and the 
less reactive identical fermions will have typical lifetimes in 
the 10 ms to 100 ms range for T around 200 nK, depending on 
species. Two exceptions are KRb and LiNa fermions, which 
have relatively small dipole moments. A lifetime on the order 
of 1 s has been achieved for like KRb fermions in the 200nK 
ranged The identical fermion lifetime will increase as l/T 
as T decreases. 



lower vibrational levels with near-universal rate constants on 
the order of 10 ~ 10 cm 3 /s. This is suggested by theoretical cal- 
culations^^]] on collisions of alkali-metal atoms with alkali- 
me tal dim ers, and by experimental measurements on Cs with 
CsjEEl Cs and Rb with RbCsEO and Cs with LiCsJSlThe ex- 
periments are at higher T where more than the lowest partial 
wave may contribute. The RbCs work^l included a specific 
theoretical calculation of the T-dependent universal rate con- 
stants summed over partial waves that agreed with the mea- 
sured quenching rate constants. Therefore, we will calculate 
the universal rate constants for the five non-reactive species in 
the lower panel of Fig. [T] and assume that they give the order 
of magnitude of the vibrational quenching rate constant in the 
F^0 limit for alkali-metal dimer molecules in states with 
v > 1 due to a collision with an alkali-metal atom or another 
dimer molecule. Future experiments can check whether this 
assumption gives good approximate magnitudes for the actual 
quenchingrate constants. 

Figure 3 shows the universal rate constants Kq and Kf(T = 
200nK) for identical bosons and fermions respectively, as well 
as the s-wave unitarity limit 



K Qu {T) = g{n%/ii){\/k), 



(5) 



where (l/k)r represents a thermal average of l/k. The uni- 
tary limit K; U (T) gives the upper bound on the rate constant. 



2.4 Effect of an electric field 

When an electric field is turned on, the reaction rate con- 
stant for like fermions increases towards unitarity with in- 
creasing field strength. This h as been explained quite simply 
by Quememer and Bohn J 31 ! 32 ! who adapted a Langevin model 
to quantum threshold conditions. This can be done more rigor- 
ously through a quantum defect model of the threshold barrier 
penetration probability. 33 Detailed analysis of free space col- 
lisions of two J = l T, dipolar molecules have been treated 
quite well in these references and othersi^HSS Here, we dis- 
cuss approximations suitable for universal ultracold collisions 
of such dipolar molecules as the electric field F is turned on. 
As collision energy increases above the threshold regime, uni- 
versal aspects of dipolar collisions also appear. 59 Since more 
than one partial wave can contribute to the collision rate as 
F increases, we include a sum over all partial waves. This 
is important for ultracold elastic collisions, although ultracold 
reaction or loss collision rates tend to be dominated by the 
contribution from the lowest partial wave in our range of elec- 
tric field strengths. 

Since scattering flux is incoming only inside the van der 
Waals radius a for such collisions and there are no scattering 
resonances, the details of the anisotropic potential for R -C a 
are not important. The collision rates are determined by quan- 
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turn scattering by the anisotropic long-range potential, which 
needs to be fully treated in a numerical approach that couples 
different \ JM0v) rotational states through the electric field and 
the anisotropic potential. For this purpose we use a coupled 
channels expansion in a \JMQv) basis. For the present we 
consider only collisions of the lowest Stark level that corre- 
lates with the J = M = molecular state for F = 0. When R is 
very large, the Stark shift in energy levels is very small com- 
pared to the spacing 2B V J to the 7=1 rotational level. The 
dominant contributions to the potential are the direct dipole- 
dipole interaction, varying as R~ 3 , and an isotropic second- 
order dispersion interaction given by a van der Waals constant 
C(,(F) that depends on F . When R becomes small enough, at 
distances R <C a, the strong dipole-dipole interaction becomes 
larger than 2B V J, and a perturbation treatment of the disper- 
sion interaction breaks down. Effectively, the two interacting 
dipoles become more strongly coupled to one another than to 
the imposed electric field F . A coupled channels expansion is 
neded to account for this changing coupling during the course 
of the collision. 

The interaction part of the Hamiltonian for the coupled 
channels expansion reads 



L H T° K 

.7=1 



■Vdd(R) + V e l(R), 



(6) 



is the Hamiltonian of a sin- 
with J ; and d 7 denoting 
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where Hf° l = B V J 2 + F • d, 

J J J 

gle molecule defined in Sec. 
the angular momentum and the dipole moment operators of 
molecule j = 1,2 respectively. The dipole-dipole interaction 

di •d 2 -3(d 1 -e R ){e R -d 2 ) 



Vdd(R) 



#3 



(7) 



depends on the projection of the dipole operator d ; on the 
axis e R = R/R connecting the center of masses of the two 
molecules. Finally, V e i(R) describes the electronic contribu- 
tion to the d ispersion interaction between molecules, given in 
"las Vei(R) = -CfR~ 6 with Cf given by an Unsold 
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Section 

approximation. 

The most complete approach is to expand the full Hamil- 
tonian, including the kinetic energy operator, in the basis 
|/iMi0v) \J2M2Ov) \£m) of eigenstates of the symmetric top 
\JjMj0v) associated with a molecule j in level v and the 
eigenstates of the angular momentum L of the relative mo- 
tion \£m). Here, i is the partial wave quantum number and 
m denotes the projection of L on the quantization axis in 
the laboratory frame. This Hamiltonian conserves the projec- 
tion of the total angular momentum on the electric field axis 
Mtot = Mi +M2 + m, and in the absence of the electric field 
also the total angular momentum Jtot = Ji + J2 + L- The 
matrix elements of the Hj and th e Vdd terms have been exten- 



sively discussed in the literature 
give explicit formulas here. 
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and we do not need to 



Typically the number of channels that have to be included 
in the quantum dynamics in the full rotational basis is very 
large. For example, the number of eigenstates with angular 
momentum /,■ < / max = 5 and partial wave quantum numbers 
(■ < ^max = 10 with Mtot = and a symmetry of the wave 
function corresponding to bosons is approximately 5700. This 
makes the coupled channels calculations highly computation- 
ally demanding. Hence, in our study of the effects of the full 
rotational expansion we apply the adiabatic approximation, di- 
agonalizing the interaction part of the Hamiltonian at differ- 
ent values of the intermolecular distance r and electric field 
F . This is motivated by the fact that the collisions of polar 
molecules in the universal y — 1 regime exhibit no resonances, 
and in free-space can be accurately modeled in the framework 
of the adiabatic potentials. 33 Our calculations with n coupled 
adiabatic channels are designated n-channel adiabatic (rota- 
tional basis), and here we consider only n = 1 . 

When the molecules are far apart, they interact very weakly 
with one another and each molecule is described by a field- 
dressed state \g(F )) with energy E g (F) and dipole d(F) along 
the field direction, found by diagonalizing the single molecule 
Hamiltonian described in Section 
Eq. |6]) can be replaced by 



2. 1 At long range H; nt in 



V int (R,9) = - 



Ce(F) 
R 6 



-d(Fj- 



l-3cos 2 



(8) 



where 9 is the angle between the axis of the electric field and 
the intermolecular axis R, and we take the zero of energy to 
be E g (F) = for a given field F. The dispersion coefficient 
C(,{F) = Cf + C™ X {F) includes a field-dependent rotational 
contribution which can be evaluated from the second-order de- 
generate perturbation theory formula (see, e.g., Refill): 



Ct{F) 



y (gim e \V dd | ei'm't) (et'm' t \ V dd \glm e ) 



Here, g and e respectively represent field-dependent ground 
and excited product states of the two molecules from the spec- 
trum of solutions to Hf°. Symmetrizing appropriately, we de- 
compose the single-molecule field dressed states into a |/M0v) 
basis, and then calculate the matrix elements of Vdd using the 
approach of Chapter 2 of Ref.^ We evaluate the sum over ro- 
tational states J and partial waves I, and sum over projections 
such that Mtot is conserved. We neglect the other degrees of 
freedom, which contribute to the non-rotational part of Cs, and 
have been calculated for some species through full electronic 
structure calculations. 47 Figure 4 gives an example of Cg ot (F) 
for the v = 0,7 = ground state of 87 Rb 133 Cs. The Figure 
shows that CJj ot (F) decreases appreciably with increasing F, 
due mostly to increasing size of the energy denominators in 
Eq. ^. 
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Fig. 4 Field-dependent rotational contribution C 1 ^ t (F) to the van 
der Waals coefficient for the interaction of two ground state RbCs 
molecules (solid blue line) as a function of electric field F, Also 
shown is the dipole moment (dashed green line). 

A description in terms of the dispersion potential in |9]) 
breaks down at small distances, when the dipole-dipole inter- 
action becomes larger than 2B V J. The characteristic length 

at which this occurs is as = (d^/ByS^^W- This length is on 
the order of 0.4a to 0.6a for the molecules we consider here. 
Consequently, since ag is already effectively in the short range 
regime where all incoming flux of particles is lost, a descrip- 
tion of collisions using Eq. ([8]l should be valid. We have car- 
ried out coupled channels (CC) calculations to test this, by us- 
ing a partial wave expansion in \lm) states with the full Hamil- 
tonian including the interaction term in the form of Eq. (jHJ. 
We call this method the CC (vdW + dipole) calculation. In 
addition, we have implemented an adiabatic approximation in 
this expansion, designated either 1 -channel adiabatic (vdw + 
dipole) or 2-channel adiabatic (vdw + dipole). 

Figure [5] shows the elastic and reactive rates calculated 
with several approximations applied to free space collisions 
of bosonic and fermionic KRb prepared in the F-dependent 
state that correlates with the v = 0,/ = state as F — > 0. For 
simplicity, we assume the same reduced mass and the dipole 
moment for bosons as for fermionic 40 K 87 Rb. The approx- 
imations are: the CC (vdw + dipole), the 1- and 2-channel 
adiabatic (vdw + dipole) and the 1 -channel adiabatic (rota- 
tional basis). All approaches predict almost identical rate con- 
stants, except for the elastic rates for bosons at high electric 
fields. In this case the rates given by the CC (vdw + dipole) 
are larger than those calculated using the adiabatic approxima- 
tions. The discrepancy stems from the ^-changing transitions 
induced by the off-diagonal elements in the dipole-dipole in- 
teraction, which are absent in the 1 -channel adiabatic approxi- 
mation. At high values of the electric field the elastic rates for 
bosons contain a significant contribution from the higher par- 
tial waves. This contribution can be included by considering 
more than one adiabatic potential. This is shown in Fig.|5]for 



bosons using the 2-channel adiabatic (vdw + dipole) approxi- 
mation, which goes much of the way to eliminating the differ- 
ence between the adiabatic and coupled channels approaches. 
Interestingly we have found numerically that calculated rates 
change by an insignificant amount if we use C(,(0) instead of 
Cs(F), Presumably this is because the van der Waals disper- 
sion term has become very small in comparison to the cen- 
trifugal and dipolar terms in the range that determines the dy- 
namics. 
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Fig. 5 Field-dependent rate constants Kn(F) and lCf(F) at 
E /Icb = 200nK versus electric field F for free space collisions of 
KRb bosons (left panel) and fermions (right panel). The figure 
shows the CC (vdW + dipole) calculation (solid lines) and 
approximations based on the 1 -channel adiabatic (rotational basis) 
(points), 1-channel adiabatic (vdw + dipole) (thin solid line), and 
the 2-channel adiabatic (vdw + dipole) approximations (dot-dash 
line). 



3 Collisions in quasi-2D geometry 
3.1 Optical lattices 

The very low kinetic energy of ultracold atoms and molecules 
allows them to be controlled and manipulated by very weak 
forces. One important control is provided by optical lattices, 
by which two counter-propagating light beams in the z direc- 
tion define a standing wave light pattern that defines a series 
of lattice cells separated by distances L in the z direction on 
the order of hundreds of nm to more than a jj.m. We con- 
sider a single isolated cell with local harmonic confinement 
along z with frequency £2/ (2n) on the order of tens of kHz. If 
kgT <C hQ,, the atoms or molecules are tightly confined to the 
ground state of the cell in the z direction, while they are free 
to move in the x — y plane. Collisions under such conditions 
are called quasi-2D collisions. On the other hand, two orthog- 
onal sets of counterpropagating light beams can provide tight 
confinement in the x and y directions, with free motion in the 
z direction. This gives rise to quasi-lD collisions. The theory 
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of quasi-2D and quasi-lD collisions, and their threshold laws, 
have been worked out in some detail,^"^ i ncluding col lisions 
of ultracold molecular dipoles in quasi-2Dj 24 l 27 ^ ^^^ 

The big advantage of reduced dimensional collisions in an 
optical lattice is the extra control one gets over collision rates 
with dipolar molecules, since the dipole length can eas- 
ily exceed the confinement length «/,. Consequently, electric 
fields can produce oriented molecules and control their ap- 
proach at long range, making it either attractive or repulsive. 
In the latter case, short range reactive or loss collisions can be 
effectively turned off, thus greatly increasing the lifetime of 
the sample. 

3.2 Collision rates in reduced dimension 

Here we will give examples of the universal collision rate 
constants for quasi-2D collisions, when the molecules col- 
lide in an electric field F in the direction z of tight confine- 
ment by a ID optical lattice. The species 40 K 8 7 Rb has al- 
ready been treated in some detail for such a caseJ^SSESl an( j 
the predicted suppression of reactive collision rates by the re- 
pulsive dipole potential barrier has been verified experimen- 
tally.^ Here we examine the transition from van der Waals 
dominated to dipole-dominated collisions for the other species 
of mixed alkali -metal dimers. The van der Waals limit applies 
when a c [ <C a, whereas the dipolar limit applies when a c [ S> «/,. 
In the latter case, collision rates approach a universal dipolar 



limit, as described in the literature ) l 70 ^ 

We use here the form of the elastic and inelastic collision 
rates from Micheli et al^\, since this makes it easy to re- 
late collisional properties in different dimensions N, where 
N = 1,2,3 refers to quasi-lD, quasi-2D, and free space col- 
lisions respectively .1^1 The "core" of the collision at short dis- 
tances R <C a is assumed to occur in full 3D geometry and to 
be represented by a normal spherical harmonic expansion in 
partial waves. The long-range part of the collision occurs in 
reduced dimensionality due to the tight confinment in 1 or 2 
dimensions with characteristic length a/, ^> a (see Fig. [3}. The 
kinetic energy E = h 2 K 2 /(2n) of free motion in the respec- 
tive (z), (x,y), or (x,y,z) directions for N = 1,2,3 is assumed 
to be <C TiQ, so the colliding molecules are confined to the 
ground state of the harmonic trap. This assumption can be 
easily relaxed to permit collisions of molecules in other trap 
levels.™ 

The wave function can be expanded in a set of partial waves 
j suitable for each dimension.^ The contributions to the elas- 
tic and inelastic rate constants for collisions with relative mo- 
mentum h K in dimension N are^l 

2 



Kf(K) = 
Kf(K) = 



nft 
M 

ith 



Sn 



gN- 



K N-2 

i-l WI 

r N-2 



(10) 



(11) 



where the factor = I /it , 2/jt , 2 for molecules colliding in 
identical spin states in N = 1,2,3 dimensions, respectively. 
The indices of the lowest partial wave are j = for bosons 
and j =1 for fermions. In 3D, j = 1 has 3 components of its 
projection M that have to be summed over. In quasi-2D with F 
along the confined direction z, j refers to the projection M of 
relative rotational angular momentum along z, with two com- 
ponents + 1 and — 1 that have to be summed over for 7 = 1; 
in quasi- ID, it refers to the symmetric (j = 0) or antisymmet- 
ric (j = 1) state propagating along z- In the case of unlike 
species both j = and 1 contribute. The contribution from the 
lowest index is dominant for ultracold universal reactive colli- 
sions in an electric field, and any additional partial waves can 
be summed over in other cases to get the total rate constants 
K ls or K ei . The loss rate for a gas with Af -dimensional density 
n is h = —K ls n 2 , where n has units of cm" 1 , cm~ 2 , and cm~ 3 
for N = 1,2,3 respectively. 

Micheli et aP^ worked out analytic expressions for any N 
for the universal K h or K ei for a van der Waals potential. The 
loss rate constants for N = 2 are Kq =2(h/n)P forM = 
bosons and Kf — 4(h/ji)Pi for fermions summed over both 
\M\ = 1 components. The prefactor 2(h/fi) or A(Ti/n) gives 
the unitarity upper bound, and the transmission function Pj = 
1 — \Sjj\ 2 gi yes me dynamical probability of getting from the 
asymptotically prepared entrance channel to small R, where 
loss occurs: 



Po 



4^n—f (K) 

Oh 



Pi 



■ (Ka) 2 , 

a h 



(12) 



where /q(k) is a complicated function that has a logarithmic 
singularity as K — > 0. However, this singularity is not impor- 
tant in the experimental nK range (it is only dominant at much 
lower energies for typical trapping geometries). For practical 
purposes /o(fc) can be taken to be on the order of unity and 
nearly independent of K for likely current experiments. 

We calculate K h here for quasi-2D geometry as a func- 
tion of electric field F using the numerical method described 
by Micheli et al. 2i similar to the method of Quemener and 
Bohn.EHlll in our case of universal collisions, we impose 
incoming-only boundary conditions on the wave function for 
R <C a, which corresponds to loss of all scattering flux that 
reaches this region from long range .^Ml The calculation thus 
requires only the van der Waals, dipolar, trapping, and cen- 
trifugal potentials, with a spherical basis expanded in spheri- 
cal harmonics at short range, switching to a cylindrical basis 
of confined states at long range. The details of short range 
"chemistry" are irrelevant, since only loss from the entrance 
channel occurs there, due to reaction of vibrational quenching. 
Since there is no back reflection to long range in the entrance 
channel, scattering resonances do not exist in this universal 
case. 

Figure |6] shows our calculations of as the electric field 
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Fig. 6 Universal loss rate constant for quasi-2D collisions at 
E/k B = 200 nK in an Q, = 2;r(50kHz) trap for bosonic (J = 0) 
6 Li 40 K, 39 K 133 Cs and 87 Rb 133 Cs (dashed lines) and K 1 * for 
fermionic 6 Li 39 K, 40 K 87 Rb and 40 K 133 Cs (solid lines). In the case 
of the nonreactive species KCs and RbCs, these rate constants are 
expected to apply to vibrational quenching for low vibrational levels 
with v > 1. The black arrow shows the unitarity limit 2fi//i for 
bosonic 6 Li 40 K. The solid points show the fit function in Eq.[T3| 
when aj > a/,. 



F is increased for a variety of species at a collision energy of 
E/kB = 200 nK for confinement in an Q. = 27r(50kHz) trap. 
These go to the universal van der Waals limits as F — >• 0. As 
we found for free space collisions, our numerical calculations 
change by an insignificant amount if we use Q (0) instead of 
C(,(F) as F increases. Consequently the van der Waals po- 
tential becomes progressively less important as F increases 
and the dipole potential becomes dominant on distance scales 
larger than the confinement length a/,. We also find that for the 
low energies and moderate dipole strengths that we explore in 
quasi-2D geometry, the lowest partial wave is the dominant 
contribution to the universal loss rate constant, so little change 
would be seen in K? in Fig. 6 if the contributions from higher 
M values were included. Elastic collisions, however, will be 
affected more, but they can be treated by simple approxima- 
tionsPm 

The 40 K 87 Rb case in Fig. [6] is the one reporte d and ex- 
plained by Micheli et a/.EHand Quemener and BohnJ^EIlThe 
other cases are new. A useful figure of merit is a rate constant 
of 10~ 7 cm 2 /s, which corresponds to a Is lifetime of a 2D 



10~ 7 cm 2 /s when d(F) is large enough. Consequently, we 
predict that stabilization of the heavier non-reactive species 
RbCs and KCs in their low vibrational levels should be possi- 
ble at dipole strength of only 1 D or less, at electric field be- 
low 10 kV/cm (see Fig. [2]). Of course, RbCs and KCs in v=0 
is nonreactive. Even a light reactive species LiK with weaker 
confinement should be stabilized in such a trap at a field on 
the order of 20 kV/cm. 

Figure [6] also demonstrates the difference between the 
bosonic and fermionic form of the same species LiK and KCs. 
As already shown for KRb, the fermionic rate is suppressed 
at zero field relative to the bosonic rate, due to the centrifu- 
gal barrier to the p-wave. At 200nK the suppression ratio 
Kf/K l Q ks 10 and 100 for KCs and LiK respectively, similar to 
the ratios in 3D (see Fig . [3j . As field increases, the fermionic 
and bosonic rate constants for the same species move away 
from the van der Waals limit and are both suppressed, even- 
tually becoming nearly equal as d(F) continues to increase. 
The figure also shows a comparison to an approximate formula 
which applies to bosons and fermions in the strong dipole limit 
where ad ^> o&, 



Kl=2-2U{Ka h ) A e- imu ^l^ 1 '\ 



(13) 



where the factor 288 comes from fitting our numerical calcu- 
lations. The factor 2 8 / 5 ps 3.0314 in the exponent is replaced 
by 2 if a ( [ and a/, are defined with the total molecular mass 
instead of the reduced mass. The agreement with the cal- 
culations shows that this expression accurately captures the 
scaling with mass and d{F) for a wide range of species. It 
seems especially good for the bosons even when a c / is only 
slightly larger than a/,. The quantity to the right of Ih/fJ, in 
Eq. ( 13 1 represents the transmission function Pq(Q.,F,E) and 
is similar to expressions for the semiclassical tunneling prob- 
ability through the long range potential barrier proportional to 
exp (— C (aj/q/; ) 275 ), wh ere C is a constant, as described by 
several authors .ISHZHsESEH Since the rate constants Kq and 
K\ s are nearly equal in the strong dipole limit, this implies 



comparable to the number 
40 K 87 RbP Smaller rate con- 



gas with density n = 10 7 cm~ 2 
achieved experimentally with 

stants correspond to lifetimes of longer than 1 s at this density. 
As the dipole strength d(F) increases, decreases below tively modest electric fields. The formula in Eq. ( 13 i is again 



Pi fa Po/2. Our fit for P (Q.,F,E) in Eq. ([13) is based on the 
calculated universal quantum dynamics and does not rely on 
any semiclassical or adiabatic approximations. It gives realis- 
tic scaling with the experimentally variable parameters: trap 
strength i2, electric field F and temperature, which respec- 
tively control confinement a/,, dipole length and the range 
of collision energy E. 

Figure [7] shows our calculations for bosonic and fermionic 
isotopes of the reactive species LiCs, which has the largest 
dipole moment of any of the mixed alkali species. Since the 
dipole moment varies slowly with vibrational quantum num- 
ber, this rate should apply to a range of low vibrational levels. 
Stabilization in quasi-2D geometry should be possible at rela- 
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Fig. 7 Universal loss rate constant for quasi-2D collisions of 
bosonic 7 Li 133 Cs (dashed lines) and fermionic 6 Li 133 Cs (solid lines) 
at E/kB = 200 nK for a range of trap strengths Q./(2n). The solid 
points show the fit in Eq. 13 when aj > a^. The arrow indicates 



the unitarity limit for bosons, 2h/fl. 



very good for this species, verifying that in the strong dipole 
limit the probability Pj decreases with increasing Q.. This is 



unlike the van der Waals limit, where Pj from Eq. ( 12 1 scales 
as Consequently Fig.|7]shows the switching between the 
two limits for dipole strengths where ad ~ ah- Again, we find 
that the universal loss rates for bosons and fermions are nearly 
the same in the strong dipole limit. 

Figure |8] shows that the A"j s in the strong dipole limit scales 
as E 2 for j = and 1 in the interesting experimental range 
of 50 nK to 500 nK, consistent with the E 2 scaling found by 
Ticknor. 12 Some departure from this scaling is evident below 
50 nK and above 500 nK in Fig. [8] In contrast, Eq. ( 12 1 shows 
that in the van der Waals limit as F — > K\ s scales as E and 
Kq is nearly independent of E. 

We have concentrated here on reactive or loss collisions in 
quasi-2D geometry, instead of elastic collisions. We have also 
calculated the latter as a function of field. Micheli et al 2& 
found that a simple unitarized Born approximation was ac- 
curate for the two lowest partial waves for ultracold KRb. As 
dipole increases it is necessary to include more than the low- 
est partial wave. Simple universal expressions for the strong 
dipole limit can be worked outJ 70 l 71 ^ The only important aspect 
of quasi-2D elastic scattering to note here is that the elastic 
rate constant continues to increase towards the unitarity limit 
(of several partial waves) as electric field increases, and thus 
the ratio of elastic to loss collison rates becomes very large 
with increasing field, due to the strong suppression of the lat- 
ter. Thus, it is hoped that fast elastic collisions can thermalize 



Fig. 8 Universal loss rate constant Kq (points) versus collision 
energy E /kg for bosonic 7 Li 133 Cs at a dipole strength of d = 0.4 au 
(1 D) for a/ (2m) = 30 kHz, 50 kHz, and 100 kHz. The solid line 
shows the fit function in Eq. d 1 3b- 



a quasi-2D gas of molecular dipoles, which are stable with 
long lifetimes relative to loss on the time scale of an experi- 
ment. 

4 Conclusion 

We have characterized the universal reactive and inelastic re- 
laxation rate constants in free space and in quasi-2D planar 
geometry of the ten mixed alkali-metal dimer polar molecu- 
lar species in the near-threshold limit of ultracold collisions 
where the lowest partial waves contribute. We consider col- 
lisions between molecules in a single state of vibration, rota- 
tion (J = 0) and internal spin, where the bosonic or fermionic 
character of the molecule leads to different rate constants. The 
universal rate constants are independent of the details of short 
range chemical interactions and depends only on the long 
range approach of the colliding species, which is subject to 
experimental control by electric and electromagnetic forces. 
The universal rates apply as long as short range chemical dy- 
namics has unit probability of chemical reaction or inelastic 
quenching that results in loss of the initially prepared ultracold 
molecules. The species LiNa, LiK, LiRb, LiCs, and KRb have 
such unit short range chemical reaction probabilities when in 
their vibrational ground state v = 0, and are expected to have 
universal rate constants in all v. The species NaK, NaRb, 
NaCs, KCs and RbCs are non-reactive in v = 0, but may have 
universal vibrational quenching rate constants when in vibra- 
tionally excited states with v > 1 . The non-universal collisions 
of the v = level of these non-reactive species are predicted to 
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have electric field-dependent scattering resonances,! 33 ! 56 ! 70 ^ 
which remain to be more fully explored theoretically and ex- 
perimentally. 

In the absence of an electric field, the universal van der 
Waals rate constants are determined only by the van der Waals 
C(, coefficient, the reduced mass of the pair, and the collision 
energy according to the known threshold laws. The impor- 
tant parameter is the van der Waals length a, which scales 
as (jUCg) 1 / 4 and thus is relatively insensitive to C(>. We use 
simple estimates of based on the electronic and rotational 
properties of the molecules at their equilibrium internuclear 
separation R e that should be adequate approximations for low 
values of vibrational quantum number v. Although the nor- 
mal second-order expansion of the dispersion energy breaks 
down when two molecules approach one another when the 
dipolar interaction exceeds the rotational spacing, the distance 
at which this happens is around half of a, where it does not sig- 
nificantly affect the reaction rate due to quantum scattering at 
distances on the order of a or larger. We find that the ultracold 
universal reactive or inelastic loss rate constants for bosonic 
mixed alkali-metal dimers should be within a few factors of 
the unitarity upper bound to the rate constant, whereas rate 
constants for fermionic species tend to be significantly smaller 
by factors of 10 to 100 around 200 nK, scaling linearly with 
temperature. 

When an electric field F is present, the collision rates vary 
with F and switch from the van der Waals universal limit to 
a universal strong dipolar limit. We show that the universal 
ultracold collision rates are not significantly affected by the F- 
dependent changes in the second-order dispersion interaction, 
but can be calculated in 3D from rotationally adiabatic po- 
tentials. We consider in detail the F-dependent universal loss 
rate constants in quasi-2D geometry when the molecules are 
oriented in the z-direction by an electric field and tightly con- 
fined in that direction by a harmonic trap such that Ka^ -C 1, 
where TiK is the collisional momentum and a/, is the harmonic 
confinement length. Our fully numerical calculations connect 
the wave function in the universal 3D short range core of the 
collision with incoming boundary conditions to the reduced 
dimensional quasi-2D wave function at long range. As in 3D 
the collision switches from the van der Waals limit to a strong 
dipole limit as F increases, where the rate constants for bosons 
and fermions of the same species are the same. We obtain a 
simple scaling formula that provides a good approximation of 
the universal rate constants for all ten species in the ultracold 
domain below 1 jj,K and for trapping frequencies in the tens 
of kHz range that are experimentally feasible. The univer- 
sal strong dipolar limit of the rate constant has simple scal- 
ing with the mass, trapping frequency, and molecular dipole 
strength. In this geometry the dipoles have long range repul- 
sive interactions, and universal collision rate constants can be 
decreased to be several orders of magnitude less than the uni- 



tarity limit, so that gases of oriented dipoles of even highly 
reactive species can be stabilized against loss collisions un- 
der practical experimental conditions that should be achiev- 
able with current techniques for cooling, trapping, and opti- 
cal lattice design and with relatively modest electric fields on 
the order of 10 kV/cm . We expect these considerations will 
also be valid for quasi- ID geometry, or collisions in "tubes" 
with tight confinement in the two directions orthogonal to the 
"tube." 

While we have considered molecular collisions for states 
with low vibrational quantum numbers v, the van der Waals 
C(, coefficient and, perhaps more importantly, the molecular 
dipole moment will vary with v. Consequently, there is some 
degree of control possible by "vibrational tuning" of these pa- 
rameters. Eventually, at very high v near the dissociation limit, 
the expectation value of the dipole moment will become very 
small, and electric field control will be lost. But there is a 
wide range of vibrational levels that should be accessible for 
experimental control, so that the v = level need not be the 
target level of STIRAP, especially if higher v levels are eas- 
ier to reach. Since the short range collision loss probability is 
maximal for universal collisions, it is not necessary to have a 
non-reactive v = species to work with, if collisions can be 
shielded in quasi-2D geometry. One could envision achieving 
a lattice of "preformed pairs" of unassociated atoms in an ar- 
ray of fully confining 3D trapping cells (zero-dimensional col- 
lisions of atoms bound to the cell), as has been proposed.! 37 ! 74 ! 
Converting such atom pairs to Feshbach molecules by magne- 
toassociation and moving the latter to a desired v by STIRAP 
should allow the study of controlled collisions upon subse- 
quent removal of the fully confining lattice by turning off some 
or all of the lattice lasers, resulting in lattices with quasi- ID 
or quasi-2D geometry or even a free space gas. A broad class 
of such collisions should have rate constants approximated by 
the universal ones described here. 
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